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Microbiotas from UC patients display 
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Microbial Ecology and Technology (LabMET), Ghent University, Coupure Links 653, B-9000 Gent, Belgium, 3 Department of 
Gastroenterology, Herlev Hospital, DK-2730 Herlev, Denmark. 

We compared fecal microbial communities derived either from Ulcerative Colitis (UC) patients in 
remission (n=4) or in relapse (n=4), or from healthy subjects (n=4). These communities were used for 
inoculation of a dynamic in vitro gut model, which contained integrated mucin-covered microcosms. We 
found that the microbiota of the 'mucus' largely differed from that of the 'lumen'. This was partly due to 
decreased mucus-associated populations of lactic acid producing bacterial populations (LAB), as LAB 
originating from UC patients had a significantly decreased capacity to colonize the mucin-covered 
microcosms as compared to those originating from healthy subjects. We found significant differences 
between the metabolomes of UC patients in relapse and remission, respectively, while the metabolome of 
patients in remission resembled that of healthy subjects. These novel findings constitute an important 
contribution to the understanding of the complex etiology of UC. 

The mucus layer lining the epithelium of the gastrointestinal tract is important for the protection of the 
intestinal epithelium in humans. Commensal bacteria have been found to colonize the colonic mucus layer, 
and previous studies have shown that the microbial community found in the colonic mucus differs from that 
of the luminal community 1,2 . Several microbial characteristics have contributed to the evolvement of the specif- 
ically selected mucosal community, including the ability of the bacteria to utilize mucin glycans as energy source 
as well as resistance to nonspecific antimicrobial peptides and specific antimicrobial immunoglobulins produced 
by the host 3,4 . Additionally, many adhesion molecules expressed by colonic bacteria have mucin glycans as 
specific epitopes 5,6 , and it has been suggested that the glycosylation pattern in mucin, hence the attachment site 
and energy source for the colonic bacteria, is an important factor for host selection of a specific mucosal 
community 3 . Lack or defects in the mucosal barrier may allow bacteria to reach the epithelium and trigger colonic 
inflammation. 

Ulcerative colitis (UC) is an idiopathic inflammatory bowel disease characterized by chronic relapsing inflam- 
mation of the colonic mucosa 7 . The etiology of UC remains an enigma, and no known infectious agent has been 
demonstrated 8,9 . It has been speculated that UC originates from a dysregulated immune response to the com- 
mensal intestinal microbiota in genetically susceptible individuals 10,11 . Human studies have revealed that UC 
patients have a colonic mucus layer that has an altered O-glycan profile and is significantly thinner than that of 
healthy subjects, which may select for a different mucosal microbial profile 12,13 . Consistently, several studies have 
shown that patients with UC have an altered bacterial microbiota 1418 . Thus, the bacterial and/or host-bacterial 
interactions may play a role in the pathogenesis of UC. 

In vitro models are well-suited to screen the adhering potency of intestinal microbes. They include adhesion 
assays to various components of the intestinal surface: e.g. intestinal mucus 1 ', mucins 20 , colonic tissue 21 or cell 
lines 22 . A drawback to such models is that they often provide only short-term information based on axenic 
cultures and thus ignore the interactions between and within the luminal and mucosal microbial communities. 
Therefore, a dynamic in vitro gut model has been developed, which simulates both the luminal and mucosal 
environment 23 . This model, named the M-SHIME was adapted from the validated Simulator of the Human 
Intestinal Microbial Ecosystem (SHIME) 24 . Compared to earlier in vitro models, the M-SHIME allows a more 
representative colonization of specific Lactobacillus sp 23 . Furthermore, high-resolution phylogenetic microbiota 
profiling shows that the simulated mucosal microbiota is, in correspondence with in vivo studies, enriched with 
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Firmicutes belonging to the Clostridial clusters IV and XlVa 24 . 
Moreover, the in vitro mucosal environment is necessary to avoid 
the wash-out of specific surface-associated microbes, which occurs in 
conventional in vitro models. Hence, the M-SHIME allows studies of 
the mucosal microbiota and the interaction between luminal and 
mucosal microbial communities. 

The aims of this study were (1) to investigate the ability of fecal 
microbiota from healthy subjects and UC patients in either remission 
or relapse to colonize the artificial mucus layer of the M-SHIME and 
(2) to elucidate microbial activity by comparison of metabolic pro- 
files of the luminal and mucosal microbial communities derived 
from UC patients and healthy subjects. Denaturing Gradient Gel 
Electrophoresis (DGGE) and quantitative Real-Time PCR (qPCR), 
both of which are culture-independent methods, were applied for 
microbiota analysis, while Liquid Chromatography - Mass Spectro- 
scopy (LC-MS) was used to analyze the potential role of extracellular 
metabolites produced by microorganisms in the lumen and mucus of 
the M-SHIME. 

Results 

Microbial activity in terms of SCFA production. SCFA were 
analyzed in the luminal content of the colonic M-SHIME vessels as 
a measure of the metabolic activity of the microbiota derived from 
healthy subjects or UC patients, either in remission or relapse 
(Table 1). Forty-two hours after inoculation with fecal samples, 
there was a significant inter- individual variability within each 
group regarding the conversion of the provided nutrients to SCFA 
(P<0.05). As a result of this variability, no significant differences 
were detectable between the three groups. However, there was a 
trend of lower acetate and higher branched SCFA concentrations 
for UC patients as opposed to healthy subjects (P = 0.138 and 
P = 0.210, respectively). Additionally, we found it noteworthy that 
the levels of carproate found in samples containing microbiota from 
patients in remission was 36 fold higher than in samples with 
microbiota from patients in relapse (P = 0.229). 

Microbial community analysis using DGGE. Comparison of 
DGGE profiles containing 16S ribosomal genes amplified from 
luminal and mucosal samples of healthy subjects and UC patients 
after 42 h colonization revealed a distinct difference between the 
dominant bacterial members of the luminal and mucosal environ- 
ment (Figure 1). The dendrogram from the Dice cluster analysis 
showed three clusters with all mucosal samples in cluster I (54.61% 
similarity), five luminal samples in cluster II (53.89% similarity), and 
seven luminal samples in cluster III (41.15% similarity). Clustering of 
the distribution of these dominant microbial species did not correlate 
to the health status of the human subjects (healthy, UC in remission 
and UC in relapse). 

Microbial community analysis by qPCR. For microbiotas derived 
from UC patients in relapse or remission, luminal and mucosal 
samples were separated from one another, while this separation 



was not equally clear for samples from healthy subjects (Figure 2, 
score plot). Roseburia, Faecalibacterium prautznitzii, and Closteri- 
diaceae/Eubacterium represented samples from the mucosal envi- 
ronment, while and B. bifidum, B. adolescentis, Bifidobacterium 
spp., Lactobacillus spp., Akk. muciniphila and Actinobacteria repre- 
sented the luminal content samples (Figure 2, loading plot). 

It was noteworthy that the preference of specific bacterial groups 
to colonize the mucosal and/or luminal compartment was clearly 
affected by the origin of the microbiota (healthy subjects, UC patients 
in remission or UC patients in relapse) (Table 2). No significant 
differences between the luminal and mucosal samples with respect 
to proportions of specific bacterial taxa were measured when fecal 
communities were derived from healthy subjects. However, in the 
vessels with communities derived from UC patients in relapse, 
mucus was colonized by significantly lower proportions of bifido- 
bacteria, B. bifidum, lactobacilli, C. coccoides group, C. leptum sub- 
group and Alistipes spp.(P = 0.01, P = 0.02, P=0.03, P = 0.05, P = 0.02, 
and P = 0.05, respectively) than found in lumen. Additionally, the 
densities of bifidobacteria, B. adolescentis, B. pseudocatenulatum, 
lactobacilli, C. leptum subgroup, Faecalibacterium prausnitzii and 
Actinobacteria derived from UC patients in remission were signifi- 
cantly lower in mucus than in lumen (P = 0.01, P = 0.001, P = 0.01, 
P = 0.007 and P = 0.05, P = 0.03 and P = 0.001, respectively). Finally, 
significantly higher mucosal counts of Roseburia spp. derived from 
UC patients in remission were measured (P = 0.05). 

Comparison of the proportions of specific taxa present in the 
luminal and mucosal compartments, respectively, from the three 
types of microbiota (healthy, UC remission, UC relapse) revealed 
that relative quantities of Bifidobacterium spp. and Lactobacillus 
spp. in mucus was significantly lower (P = 0.05 and P = 0.002, respect- 
ively) in communities derived from UC patients in relapse than in 
those derived from healthy subjects, whereas the relative quantity of 
Clostridiaceae/Eubacterium was significantly higher in the lumen 
(P = 0.04) in communities from UC patients in relapse than in those 
from healthy subjects (Table 3). 

The bacterial growth rate of selected bacterial taxa in the lumen 
was calculated from the slope of the exponential phase of the growth 
curves. The growth rate was increased during the 42 h incubation in 
the lumen for all the examined bacterial taxa. No significant differ- 
ence in the growth rate for each bacterial taxon was found when 
comparing the three groups (healthy subjects, UC in remission or 
UC in relapse) (Table S3, supplementary data). 

Metabolite detection and separation. PCA of the metabolites as 
detected by LCMS revealed a difference between samples taken 
from healthy subjects and UC patients in relapse, respectively. The 
grouping was present in samples from mucus (Figure 3, score plot) as 
well as from lumen (Figure 4, score plot). No difference was observed 
between UC patients in remission and healthy subjects (Figure SI 
and S2, supplementary data). However, a clear separation between 
score plots for UC patients in relapse and remission, respectively, was 
seen (Figure S3 and S4, supplementary data). 

Based on the PCA loading plots (Figure 3 and 4), metabolites, 
which were tentatively causative for the difference observed between 
UC patients in relapse and healthy subjects, were identified (Table 4 
and 5). Loadings, which in extracted ion chromatograms showed 
clear chromatogram peaks (data not shown) that were present in 
significantly different levels (P<0.05) in the two groups were further 
investigated. Substances in the Human Metabolome Database, which 
corresponded to the found mass/charge ratio (M LCMS ) were iden- 
tified. However, a few of the given mass/charge values had more than 
one possible match (Table 4 and 5). 

Discussion 

In this study, we applied a recently developed dynamic in vitro 
gut model, the M-SHIME, to investigate differences between the 



Table 1 The average ( ± SEM) 


absolute SCFA levels (mM) in the 


luminal content of the M-SHIME 


units, 42 h after 


inoculation with 


fecal samples of different human subjects: healthy, UC remission 


and UC relapse 


n = 4) 








Healthy 


UC remission 


UC relapse 


Acetate 


54.0 ± 3.8 


40.7 ± 4.9 


46.1 ± 6.0 


Propionate 


7.3 ± 0.3 


6.2 ± 3.2 


9.4 ± 1 .4 


Butyrate 


23.6 ± 5.3 


18.9 ±4.2 


23.8 ± 3.4 


Valerate 


0.6 ± 0.6 


1.6 ± 0.8 


0.7 ± 0.6 


Caproate 


1.5 ± 1.5 


3.6 ± 2.8 


0.1 ±0.1 


Branched SCFA 


2.0 ± 1.6 


5.1 ± 1.2 


3.2 ± 1.3 


Total SCFA 


89.1 ± 2.5 


76.2 ± 4.6 


83.3 ± 7.6 
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Figure 1 | Dice cluster analysis of universal DGGE gel profiles from in vitro luminal and mucosal communities of the M-SHIME colonized with 
samples derived from healthy subjects,or from UC patients in remission or relapse. The luminal samples are indicated by stars (*) and the mucosal 
samples are indicated by full circle (•). The dendrogram can be divided into three clusters: Cluster I: mucosal samples (54.61% similarity). Cluster II: 
luminal samples (53.89% similarity). Cluster III: luminal samples (41.15% similarity). Metric scale indicates degree of similarity in percentages. 
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Figure 2 | Principal Component Analysis of the quantitative PCR measurements illustrated by. PCI and PC2 (29.9% and 21.3% of explained variance, 
respectively). Score plot showing the M-SHIME luminal (A) and mucosal ( O ) communities. Sources of the communities are indicated by green for 
healthy subjects, blue for UC patients in remission and red for UC patients in relapse. Loading plot indicating each of the measured bacterial taxa as 
determined by quantitative Real-Time PCR. 1. B. bifidum; 2. B. adolescentis; 3. B. pseudocatenulatum; 4. Bifidobacterium spp.; 5. Lactobacillus spp.; 6. C. 
Upturn subgroup; 7. C. coccoides group; 8. F. prausnitzii; 9. Dcsulfovibrio spp.; 10. Akk. muciniphila; 1 1. Firmicutes; 12. Bacteroidetes; 13. Roseburia spp.; 14. 
Bacteroides spp.; 15. Alistipes spp.; 16. Actinobacteria; 17. Bac. fragilis group; 18. Clostridiaceae/Eubacterium. 
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Table 2 | Preference of bacterial taxa to colonize the mucosal compartment, expressed as ratio of the relative quantities in mucosal and 
luminal compartments of the M-SHIME 



Ratios mucus/lumen (%)" 



Bacterial taxQ 


Health 


subjscts 


UC remission 


UC relapse 


Firmicutes 


100.36 


(±1.23) 


1 00.46 


>— ' ■°4) 


99.85 


[±0.89) 


• Clostridiaceaec b /Eubacterium 


103.18 


;±29.64) 


1 1 ~7 A ~7 
1 1 / .4/ 


± 1 V.4VJ 


"70 

/V .Do 


l-t-"7 OOt 


• Clostridium leptum subgroup 


96.33 


;±8.43) 


83.81 


[±5.1 2)* 


82.39 


[±3.35)* 


° Faecalibacterium prausnitzii 


102.62 


±13.35) 


83.75 


[±4.22]* 


93.91 


1+4.97) 


• Clostridium coccoides group 


88.47 


;± 10.94) 


74.10 


[± 1 1 .74) 


66.81 


[±10.73)* 


° Roseburia spp. 


108.01 


±7.70) 


137.83 


± 1 2.28)* 


146.06 


[±41.24] 


• Lactobacillus spp. 


94.70 


;±8.85) 


81 .32 


[±2.77)** 


68.59 


[±7.64)* 


Bacteroidetes 


99.98 


±2.84) 


106.52 


[±9.69) 


94.36 


[±7.57] 


• Bacteroides spp. 


98.44 


;±2.87) 


1 12.93 


±9.33) 


85.28 


[±5.06) 


° Bac. fragilis group 


93.95 


±7.55) 


107.24 


±14.33) 


75.48 


[±8.88) 


• Alistipes spp. 


65.05 


;i7.58) 


68.12 


;± n.76] 


71.91 


(±8.47)* 


Actinobacteria 


76.29 


13.79) 


55.39 


;2.88)*** 


52.01 


[16.96] 


• Bifidobacterium spp. 


94.60 


;±3.12) 


91.24 


±1.53)** 


84.63 


[+4,28)** 


° 8. bifidum 


93.99 


±3.47) 


106.37 


±6.97) 


84.48 


[±3.35]* 


° 8. adolescentis 


89.52 


±4.61) 


83.68 




67.61 


[±16.58] 


° 8. pseudocatenulatum 


101.27 


;±4.68) 


79.79 


±3.55)** 


94.60 


(±15.00) 


Proteobacteria 














• Desulfovibrio spp. 


101.19 


(±1.75) 


101.22 


(±4.53) 


92.74 


(±3.89) 


Verrucom icrobia 














• Akkermansia muciniphila 


85.48 


(±6.90) 


73.15 


(±12.19) 


86.18 


(±5.10) 



"Ratios (%) calculated as 1 00*mucosal samples/luminal samples, hence values below 1 00% indicates low adherence capacity to mucus, while values above 1 00% indicate high adherence capacity to 
mucus. 

Includes Clusters I, III, IV, XlVa, XlVb. 

Asterisks (*) indicate significant differences between mucosal and luminal samples (*P<0.05 and **P<0.01). All calculated data are means ± SEM. 



intestinal microbial ecosystems of healthy subjects and UC patients, 
either in relapse or in remission. This model allows investigation of 
differences within the luminal content as well as at the artificial 
intestinal mucosal surface 23,24 . The impact of the human host on 
the microbial composition is eliminated, thus allowing focus on 
intrinsic features of the gut microbial populations. However, as 
any model, the M-SHIME has its limitations, and it should be noted 



that the bacterial communities colonizing an in vitro fermentor can 
never be made completely identical to those present in humans due 
to e.g. possible loss of species during freezing, and altered availability 
of nutrients leading to changes in the competitive advantages of other 
species. 

Still, in line with recent in vivo 25,26 and in vitro studies 27,28 , we have 
demonstrated that also in the M-SHIME model, the in vitro mucosal 



Table 3 | Relative abudance of bacteria derived from UC patients compared to healthy subjects 

Lumen 

UC remission 



Mucus 



Bacterial taxa 



UC relapse 



UC remission 



UC relapse 



Firmicutes 


98.99 


±2.34) 


101.92 


±0.52) 


99.02 ( 


±0.96) 


101.43 ( 


±0.40) 


• Clostridiaceaec°/Eubacterium 


1 1 1 .60 


±13.44) 


142.63 


±8.74)* 


1 3 1 .70 ( 


±9.89) 


1 20.27 ( 


±13.98) 


• Clostridium leptum subgroup 


105.94 


±9.04) 


120.45 


±3.06) 


95.21 ( 


± 1 0.49) 


106.03 ( 


±4.53) 


° Faecalibacterium prausnitzii 


120.28 


±7.78) 


125.79 


±3.57) 


103.39 ( 


±2.32) 


122.12 ( 


±5.42) 


• Clostridium coccoides group 


97.68 


±14.02) 


1 19.86 


±6.22) 


69.89 ( 


±23.95) 


94.02 ( 


±16.97) 


° Roseburia spp. 


103.12 


±7.65) 


72.29 


±1 1.66) 


133.15 ( 


±19.83) 


89.96 ( 


±24.74) 


• Lactobacillus spp. 


109.74 


±3.52) 


80.06 


±12.88) 


96.27 ( 


±4.42) 


57.99 ( 


±8.93)*** 


Bacteroidetes 


95.33 


±1 1.68) 


107.10 


±4.71) 


98.92 ( 


±6.41) 


100.04 ( 


±4.25) 


• Bacteroides spp. 


93.35 


±1 1.66) 


1 10.06 


±6.83) 


1 04.27 ( 


±6.83) 


94.82 ( 


±6.63) 


° Bac. fragilis group 


94.12 


±15.12) 


1 13.49 


±8.76) 


103.68 ( 


±14.08) 


88.68 ( 


±6.82) 


• Alistipes spp. 


100.97 


± 1 1 .46) 


1 1 1.57 


±6.01) 


100.60 ( 


± 1 9.04) 


1 18.90 ( 


±18.49) 


Actinobacteria 


101.99 


±7.42) 


1 00. 1 0 


±3.65) 


74.06 ( 


±8.91) 


68.72 ( 


±23.89) 


• Bifidobacterium spp. 


99.97 


±3.72) 


97.65 


±2.65) 


96.53 ( 


±4.95) 


89.04 ( 


±4.22)* 


° 8. bifidum 


71.84 


±24.63) 


79.50 


±1 1.98) 


75.73 ( 


±24.27) 


71.53 ( 


±1 1.61) 


° 8. adolescentis 


1 12.86 


±4.80) 


89.80 


±15.22) 


106.56 ( 


±5.15) 


76.15 ( 


±23.01) 


° 8. pseudocatenulatum 


1 13.94 


±7.96) 


96.16 


±22.55) 


99.55 ( 


±8.1 1) 


90.74 ( 


±18.72) 


Proteobacteria 


















• Desulfovibrio spp. 


1 17.09 


±17.88) 


120.72 


(±5.61) 


1 17.92 ( 


±18.33) 


1 1 1 .03 ( 


±4.90) 


Verrucom icrobia 


















• Akkermansia muciniphila 


90.88 


±19.54) 


79.57 


(±17.82) 


85.1 1 ( 


±24.57) 


73.99 ( 


±14.50) 



All numbers are average ± SEM of the four samples in each UC group of either lumen or mucus. The relative abundance values (%) are normalized to healthy controls (set to 1 00). 
includes Clusters I, III, IV, XlVa, XlVb. 

Asterisks (*) designate a significant difference from the healthy control group (*P < 0.05 and ***P<0.001 ). 
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Figure 3 | Principal Component Analysis of LCMS data from the mucosal M-SHIME samples after colonization with microbiota derived from either 
healthy subjects (O) or UC patients in relapse (A). Score (left) and loading (right) plots are shown. PCI and PC2 explain 29.1% andl8.6% of the 
variance, respectively. Numbers from lto 10 in the loading plot correspond to metabolite candidates enriched in the M-SHIME mucin-coveres 
microcosms.. Masses of the candidates are given in Table 4. 



microbial community differs from the luminal one (Figure 1), with a 
tendency of specific butyrate-producing bacteria (e.g. Roseburia 
spp.) being abundant members of the mucosal microbiota (Table 2 
and Figure 2). In samples from healthy people, the difference 



between mucosal and luminal populations was more clearly seen 
by DGGE than by qPCR, probably reflecting that by qPCR only 
the subpopulations targeted by the selected primers are seen. The 
luminal and mucosal bacteria have previously been demonstrated to 



o • 



13 



12-11 



14< 



o 

o o 




• 16 



o o 



• 15 



Figure 4 | Principal Component Analysis of the LCMS data from the luminal M-SHIME samples after colonization with microbiota derived from 
healthy subjects (O) or UC patients in relapse (A). Score (left) and loading (right) plots are shown. PCI and PC2 explain 26.2% and 21.8% of the 
variance, respectively. Numbers from 1 1-1 7 in the loading plot correspond to metabolite candidates that are enriched in the luminal M-SHIME samples. 
Masses of the candidates are given in Table 5. 
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No. 


Metabolite candidate 


Mrca (Da)/RT (min 


Mlcms (Da) 


Mhmdb (Da) 


Error (mDa) 


P values 


i 
i 


Mnnni i n cfiti i rntori rfitH/ nnnc 1(1 R 1 

1 V ILJI IUUI I5UIUI vJICVlJ IUMy ULIUj I V_. 1 (Jl 


281 .5/6 


281 24R569 


ZO 1 .Z400Z/ 


U.VJOO 


r ^u.uo 


o 
z 


Polyunsaturated fatty acids (CI 8) 


279.5/6 


279.232812 


970 9T9Q71 

Z/ 7 .ZJ Z 7/ 1 


0 1 ^0 


r^u.uj 


o 
o 


Secondary bile acids 


391.5/4 


391.285401 


J7 1 .ZOJ4UU 


o on l 


r^u.uj 


4 


Salicyluric acid (drug) 


1 94.5/2 


1 94.046042 


194.045883 


0.159 


P<0.001 


5 


Monounsaturated fatty acids (CI 8) 


282.5/6 


281.248569 


281.248627 


0.047 


P<0.05 


6 


A ■ 1- 1 ■ -l/l - | 1 - ■ % r\ 

Aminosalicylic acid (drug metabolite) or 3- 


152.5/1 


152.035358 


1 59 035309 


0.049 


P<0.01 




Hydroxyanthranilic acid (oxidation product of 














tryptophan metabolism) 












/ 


Product of phenylalanine metabolism 


1 AQ "i/T 
1 7 . o / o 


1 AQ C\A(Y7hA 


1 AO MAMRHA 
1 4V.UOUOUO 


U.U4Z 


r<u.u l 


Enric 
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INO. 


Metabolite candidate 


Mrca (Da) 


Mlcms (Da) 


Mhmdb \W) 


Error (mDa) 


D 1 

P values 


8 


Tryptophan 


203.5/2 


203.082491 


203.082596 


0.105 


P<0.05 


9 


Phenylalanin 


164.5/1 


164.071619 


164.071701 


0.082 


P<0.005 


10 


Hydroxyphenyllactic acid or 3-(3-hydroxyphenyl)-3- 


181.5/2 


181.050629 


181.050629 


0.531 


P<0.05 




hydroxypropanoic acid 












Numb 


ersfrom 1-10 correspond to the metabolite marked in the loading plot Fi 


lure 3. Mpca designates mass 


faken from PCA; RT designates retention time bucket 


of the PCA; M LC ms 


designates mass taken 
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display different roles in the host, and it has been proposed that the 
mucosal microbiota is more involved in interaction with the epithe- 
lial and immune cells than the luminal microbiota, because it resides 
closer to the intestinal epithelial cells 29 ' 30 . Hence, an altered mucosal 
microbial community may play an important role in dysregulated 
immune responses. 

In spite of a limited sample size, we found a number of statistically 
significant differences between microbial communities derived from 
healthy people, UC patients in remission, and UC patients in relapse, 
respectively. For samples from UC patients a significantly lower 
proportion of lactobacilli and bifidobacteria residing in the mucus- 
microcosm than in the lumen was observed, while this was not the 
case for healthy microbiotas (Table 2). Notably, this was observed for 
UC patients in remission as well for those in relapse, indicating that 
the impaired ability of the lactic acid bacteria to adhere to mucus is 
present also when UC patients are free of severe symptoms. 
However, when directly comparing 'mucosal' populations from 
M-SHIME colonized with microbiota from UC patients to those 
colonized with healthy microbiotas, only patients in relapse had a 
significantly lower amount of bifidobacteria and lactobacilli than 



found in the 'healthy' samples (Table 3). This could be due to the 
fact that microbiota from UC patients in relapse in general contained 
lower amounts of lactic acid bacteria than microbiota derived from 
either healthy people or UC patients in remission. Species or strain 
specific mucus adhesion promoting proteins have been reported in 
several bifidobacteria and lactobacilli 6,3135 . The expression of adhe- 
sion molecules may be changed in the lactic acid producing bacteria 
derived from UC patients, hence their inability to colonize the mucus 
in vitro. However, we would expect that host-induced differences in 
bacterial gene expression profiles would no longer prevail after 
42 hours in the gut model. Alternatively, an altered composition of 
the species of lactobacilli and bifidobacteria derived from UC 
patients may explain the decreased ability of these groups to adhere 
to the mucin-microcosms of the M-SHIME. A third explanation may 
be that within a given species, the strains of lactic acid bacteria 
present in UC patients are less capable of adhesion to mucins than 
their counterparts present in healthy subjects, perhaps due to minor 
genetic differences in adhesion-relevant genes. 

Bifidobacteria and lactobacilli are believed to play important roles 
in promoting intestinal health 36 43 . The observed depletion of bifido- 
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15 Tryptophan 
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165.5/3 
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bacteria and lactobacilli in UC patients in relapse, which is in line 
with a number of previous reports 44 46 , is likely to have a consequence 
for colonic health. It is, however, not possible from the given data to 
conclude whether the impaired adhesion capacity of lactic acid bac- 
teria derived from UC patients is part of the cause or part of the effect 
of the disease. It may be that the mucus layers in the intestines of the 
UC patients have selected for populations of lactic acid bacteria, 
which are less optimized to adhere to the 'healthy' mucin of the 
M-SHIME. 

Also the ability of the butyrate-producing clostridial groups C. 
coccoides and C. leptum to colonize the mucin-covered microcosms 
of the M-SHIME was found to be lower in bacterial communities 
from UC patients than in those originating from healthy subjects, 
although for UC patients remission, this was only statistically sig- 
nificant for the C. leptum group (Table 2). This may be linked to the 
lower abundance of lactobacilli and bifidobacteria in the micro- 
cosms, since these species produce acetate and/or lactate, which is 
subsequently utilized by the given clostridial groups. In contrast to 
the C. coccoides group and the C. leptum subgroup, the buturate- 
producing Roseburia spp. in microbiotas originating from UC 
patients demonstrated a high presence mucin-microcosms 
(Table 2), indicating that these species, which are known to be able 
to degrade complex carbohydrates 47 , are not dependent on the pres- 
ence of lactate and acetate as seen for their clostridial relatives. 
Hence, the low abundance of the clostridial groups within the 
mucin-microcosms may create a free ecological niche allowing for 
excessive colonization of Roseburia spp. 

Metabolomic (LCMS) analysis revealed that the microbiotas from 
healthy subjects and from UC patients in remission displayed very 
similar metabolisms, while the metabolism of bacterial communities 
from UC patients in relapse was clearly different from these two 
groups (Figure 3, 4, SI, S2, S3 and S4). We found that metabolism 
of bile acids, tryptophan and phenylalanine were altered in luminal 
as well as mucosal samples derived from UC patients in relapse as 
compared to healthy subjects (Table 4 and 5). The presence of drugs 
and drug metabolites in samples from UC patients in relapse was 
expected, since all four patients received mesalazine ^-aminosa- 
licylic acid). 

In the human body, most bile acids are actively reabsorbed from 
the ileum and returned to the liver. However, a small fraction of bile 
acids escape enterohepatic circulation and enter the colon where 
bacteria metabolize the bile acids primarily by deconjugation and 
oxidation of hydroxyl groups. Bile salt hydrolysis is carried out by 
a broad spectrum of intestinal bacteria 48,49 . However, the specificity 
of hydroxysteroid dehydrogenase (HSDH) varies depending on bac- 
terial species and may either be specific for the 3-, 7-, and/or 12- 
hydroxy groups of bile acids leading to different secondary bile acids. 
Specific HSDHs have primarily been detected in a number of 
Clostridium and Eubacterium species belonging to either cluster I, 
IV (C. leptum subgroup), XlVa (C. coccoides group) or Xlb 49,50 . As 
we demonstrated significantly higher levels of Clostridiaceael 
Eubacterium in the luminal samples from UC patients in relapse 
than in healthy subjects (Table 3), this could partly explain the 
observed enrichment of secondary bile acids. High levels of second- 
ary bile acids such as deoxycholic acids (primarily produced by spe- 
cies belonging to Clostridiaceae) can be detrimental for colon 
health 51 . 

Phenylalanine was more abundant in mucosal samples from the 
M-SHIME colonized with microbiota from healthy subjects, than in 
those colonized with samples from UC patients (Table 4). Products 
of phenylalanine that are normally metabolized by intestinal bacteria 
undergo a variety of processes in the body, where they may be detoxi- 
fied by either glucuronide or sulphate conjugation in the gut mucosa 
and liver, or remain unabsorbed and voided in feces 52 . Production of 
phenolic compounds by proteolytic and peptidolytic activities of the 
intestinal bacteria has previously been associated with various dis- 



eases including cancer and UC 52 . Previous studies have shown that 
species belonging to the genera Clostridium and Bacteroides have a 
high capacity to ferment phenylalanine to phenolic compounds such 
as phenylpropionate, phenylacetate and/or phenyllactate 53,54 . Here, 
we observed higher levels of Clostridia in samples derived from UC 
patients than in samples from healthy subjects (Table 3). Hence, we 
speculate that the enrichment of products of phenylalanine metabol- 
ism in mucosal and luminal samples derived from UC patients in 
relapse could be a result of the higher abundance of these bacterial 
groups (Table 4 and 5). 

Also tryptophan levels were higher in samples from the M-SHIME 
colonized with microbiota from healthy subjects than when it was 
colonized with microbiota from UC patients in relapse (Table 4 and 
5). In line with this, previous studies have revealed that after fer- 
mentation using inocula from healthy subjects, no products of tryp- 
tophan metabolism could be detected in the system 52,54 . 

Finally, we observed that unsaturated fatty acids (FA) were sig- 
nificantly higher in mucus compartments colonized with bacteria 
from UC patients in relapse than in those colonized with 'healthy' 
bacterial communities (Table 4). Previous metabolomic research has 
revealed changes of FA in serum of UC patients 55 , and has suggested 
that these changes could be caused by increased endogenous biosyn- 
thesis of polyunsaturated fatty acids (PUFA) and by increased lipo- 
lysis initiated by cytokines during the inflammatory response. Based 
on the present work, we speculate that the altered levels of FA in 
serum UC patients in relapse compared to healthy subjects could be 
originating from the intestinal microbial metabolism. 

In conclusion, our data suggest that the significantly altered intest- 
inal bacterial community present in UC patients in relapse results in 
a significantly altered bacterial metabolic profile. The observed bac- 
terial alterations are suggested to result in increased metabolism of 
phenylalanine and tryptophan in microbial communities from UC 
patients in relapse. While UC patients in remission display micro- 
biotas and metabolomes very similar to those of healthy subjects, the 
lactic acid bacteria present in patients in remission are, similar to 
those from the relapse samples, significantly impaired in their ability 
to adhere to the mucus microcosms of the M-SHIME. This may be 
due to a different composition or a different genetic makeup of the 
lactic acid bacteria present in all UC patients, which may play a role 
in the etiology of this disease. We suggest that probiotic therapy for 
UC patients should not exclusively aim at increasing the amount of 
lactic acid bacteria present in the gut, which has previously proved 
helpful 56 , but also at replacing the existing lactic acid bacteria with 
other strains/isolates with better capacity for mucosal adhesion. 

Methods 

Human volunteers and clinical characteristics of UC patients. Fecal samples were 
obtained from 8 patients with UC and 4 healthy controls 57 . Within the UC group, 4 
patients were in clinical remission and 4 patients had active disease at the time of 
sampling according to clinical and endoscopical criteria 58 . The study was performed 
in accordance with the Second Helsinki Declaration, reported to the Danish Data 
Protection Agency and approved by the Regional Ethics Committee. Written 
informed consent was obtained from each participant under a protocol approved by 
the Danish National Committee on Biomedical Research Ethics. All four patients 
with active UC were treated with mesalazine (5-aminosalicylic acid). Two patients 
with inactive UC received mesalazine, one patient received olsalazine (6- 
hydroxybenzoate) -salicylic acid), and one received no antiinflammatory treatment. 
None of the participants had been treated with antibiotics for at least 2 months before 
enrolment and there was no significant difference (P — 0.32) in the mean age of the 
participants when comparing the 3 groups. 

Sample collection and processing. Stool samples were collected in airtight containers 
and stored at 4°C (limited storage time was encouraged 59 ) until delivery to the 
laboratory. Feces were homogenized in glycerol to give a 25% feces/glycerol slurry. 
This was performed in an anaerobic cabinet (Macs Work Station, Don Whitley, 
containing 10% H2, 10% C02, and 80% N2). The processed samples were stored at 
— 80 U C until further analysis. 

Growth medium and chemicals. Unless stated otherwise, chemicals were obtained 
from Sigma (Bornem, Belgium). The M-SHIME feed contained 1.0 g/1 
arabinogalactan, 2.0 g/1 pectin, 1.0 g/1 xylan, 3.0 g/1 starch, 0.4 g/1 glucose, 3.0 g/1 
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yeast extract, 1.0 g/1 peptone, 4.0 g/1 mucin, and 0.5 g/1 cystein. Pancreatic juice 
contained 12.5 g/1 NaHC03, 6.0 g/1 bile salts (Difco, Bierbeek, Belgium) and 0.9 g/1 
pancreatin. Mucin agar was prepared by boiling autoclaved distilled H 2 0 containing 
5% porcine mucin type II and 1% agar. The pH was adjusted to 6.8 with 10 M NaOH. 

M-SHIME. Analysis in the M-SHIME dynamic gut model was carried out as 
previously described 27 at two different occasions with six participants for each 42- 
hours run (two healthy, two UC patients in remission and two UC patients in relapse). 
The position (vessel) of the inocula from either healthy subjects or UC patients was 
changed for each run. 

Microbial activity in terms of short-chain fatty acids (SCFA). Acetate, propionate, 
butyrate, valerate, caproate and branched SCFA (isobutyrate, isovalerate and 
isocaproate) were measured as described previously 60 . 

Extraction of bacterial DNA. Before extraction of the mucosal samples, the samples 
were heated for 15 min at 55 'C to make the agar soluble. Subsequently, DNA was 
extracted from thawed samples using the QIAamp DNA Stool mini kit (Qiagen, 
Hilden, Germany) with a bead beater step in advance, as previously described 61 . For 
each sample, DNA was extracted in duplicates. The purified DNA was stored at 
— 20°C until use. 

PCR amplification for DGGE. Aliquots (10 uL) of purified DNA (5 ng/ul of pooled 
DNA from the duplicate DNA extractions) were applied to the following to give a 
50 uL PCR reaction mixture: 20 ~L of 5 PRIME MasterMix (2.5 X ) (VWR & Bie & 
Berntsen) and 10 pmol of each of the primers (Eurofins MWG Synthesis GmbH, 
Ebersberg, Germany). Primers HDA1-GC/HDA2 62 targeting 16S rRNA genes from 
all bacteria were used in a touchdown PCR as previously described 63 . 

Analysis of luminal and mucosal microbiota by DGGE. DGGE was carried out as 
described previously 64 using a DcodeTM Universal Mutation Detection System 
instrument and gradient former model 475 according to the manufacturer's 
instructions (Bio-Rad Labs, Hercules, California). The 9% polyamide gels were made 
with denaturing gradients ranging from 25% to 65%. The 100% denaturant solution 
contained 40% formamide and 7 M urea. Thirteen microlitres PCR products were 
mixed with 3 uL loading dye before loading. Gels were run in 1 x TAE at 60 °C for 
16 h at 36 V, 28 mA, stained with ethidium bromide for 15 min, destained for 
20 min, and viewed by UV-B trans illumination at 302 nm (Bio-Rad). The 
BioNumerics software, version 4.60 (Applied Maths, Sint- Martens -Latem, Belgium) 
was used for identification of bands and normalization of band patterns from DGGE 
gels. DGGE gels were normalized by an assigned marker (developed in our 
laboratory). A cluster analysis was performed based on Dice coefficient of similarity 
(weighted) using the unweighted pair group method and the arithmetic averages 
clustering algorithm. 

Quantitative PCR assay conditions. QPCR was performed on an ABI Prism 
7900 HT from Applied Biosystems. The amplification reactions were carried out in a 
total volume of 1 1 uL containing; 5.50 uL (SYBR* Green, Applied Biosystems), 
primers (each at 200 nM concentration) (Eurofins MWG Synthesis GmbH), 2 uX 
template DNA, and Nuclease-free water (Qiagen) purified for PCR. The amplification 
was carried out as previously described 65 . DNA (5 ng/ul) from the duplicate DNA 
extractions of each sample was used for the qPCR. 

Quantitative PCR primer and data handling. The primers specific to regions of the 
16S rRNA genes of 20 selected bacterial taxa are listed in Table SI (supplementary 
data). The taxa were selected based on an expected relevance to human health, and the 
specificity of the primers under the given conditions has recently been verified 66 . The 
relative quantities of gene targets encoding gene sequences of the bacterial taxa were 
calculated using 2 DeltaCt , assuming primer efficiency at 1.0. Delta Ct is the Ct-values of 
the bacterial target normalized against Ct-values of the total bacterial population in a 
sample. Ct is the threshold cycle calculated by the ABI software (SDS version 2.2; 
Applied Biosystems, Foster City, California, USA) as the PCR cycle, where 
amplification signal exceeds the selected threshold value, also set by the ABI software. 
Prior to the quantification, standard curves were created using serial 10-fold dilutions 
of bacterial DNA extracted from one of the M-SHIME samples for all primer sets. 
Analysis of the standard curves allowed verification of PCR efficiency for the chosen 
PCR conditions. All primers were tested to confirm sensitivity and specificity using 
DNA from pure bacterial species (Table S2, supplementary data). The detection limit 
was set to 0.001% of the quantity of the total bacteria. Bacterial targets that could not 
be detected or were below the detection limit were set to one half of the detection limit 
for further calculations. 

Metabolite detection & identification. Luminal samples from the M-SHIME were 
thawed on ice and subsequently centrifuged at 3000 g for 15 min at 4°C. Cold 
methanol (around — 80 C C) in a ratio 1 : 1 was added to the supernatant. Metabolites 
from the mucosal samples, due to the semi-solid form, were extracted with 2 ml of 
cold 50% methanol and centrifuged at 3000 g for 15 min at 4 C. This procedure was 
repeated twice. All of the samples were stored at — 80°C until further analysis and 
centrifuged at 15000 g for 5 min at 4°C just before LCMS analysis. The analysis of 
metabolites was conducted as follow: a Dionex Ultimate 3000 RS liquid 
chromatigraph (Dionex, Germering, Germany) coupled to a Bruker maXis time of 
flight mass spectrometer equipped with an electrospray interphase (Bruker Daltonics, 



Bremen, Germany). Analytes were separated on Kinetex pentafluorophenyl column 
50 x 2.10 mm, 2.6 um, 100 A (Phenomenex, USA), using the solvent system: A, 
10 mM ammonium formate pH 3.5, and B, acetonitrile. Solvent programming was 
isocratic 0% B to 0 min followed by a linear gradient to 100% B at 7 min and 100% B 
at 8 min. Flow rate was 0.25 ml/min at 0 min, increased to 0.4 ml/min at 7 min. 
Solvent composition and flow were returned to initial conditions at 8.2 min. The oven 
temperature was 40°C. Injection volumes were 1 ul. The following electrospray 
interphase settings were used: nebulizer pressure 2 bar, drying gas 10 1/min, 200 C, 
capillary voltage 4000 V. Scan range was from 50 to 800 m/z. The main focus of this 
study was acids, thus samples were analyzed in negative mode. 

To identify metabolites seen to separate given groups, the Human Metabolome 
Database (HMDB; Wishart et al. 2009) was used. Molecular formulas of the meta- 
bolites from the LCMS analysis were generated based on exact mass and isotopic 
pattern with Bruker Daltonics Software. Identification of metabolites was based on a 
very low measurement error (Tab.l, Tab. 2; external & internal calibration and lock 
mass was used), compared to the data given in HMDB and further references to the 
microbial metabolism in the GI. In addition hydrophobic properties of the molecule, 
indicated by the retention time in the chromatograph, were taken into consideration. 

Statistical analysis. Statistical analysis of the qPCR data was performed using 
OriginPro software (version 8.1; OriginLab Corporation, Northampton, USA). 
Normality and homogeneity of variances of the qPCR data were assessed using the 
Kolmogorov-Smirnov and the Levene's test, respectively. Log-transformations were 
performed for data that did not meet these criteria. T-testing was apllied to identify 
significant differences in colonization ability using the ratio (proportion) of the given 
bacterial taxon present in the two communities (lumen and mucus) for the three 
groups (healthy, UC patients in remission and UC patients in relapse). Univariate 
ANOVA was used to determine significant differences between specific bacterial taxa 
comparing the three groups either from lumen or mucus. Where ANOVA indicated a 
significant difference, Fisher's least significant different test was used. The 
nonparametric Kruskal-Wallis test was used for datasets, which did not have 
homogeneity of variance or were not normally distributed even after log- 
transformation. Univariate ANOVA was applied to confirm that there was no age 
differences between the three disease groups. Tests were considered statistically 
significant if P-values lower than 0.05 were obtained. Principal Component Analysis 
(PCA) of qPCR data was carried out using LatentiX® data analytical software (version 
2.00, The Mathworks Inc., Copenhagen, Denmark). Differences in metabolite profiles 
were evaluated by PCA using Profile Analysis 2.0 by Bruker Daltonics. Data was 
grouped into buckets of 1 min and 1 m/z differences in the range from 0.5 to 9 min 
and 50 to 800 m/z and normalized by the sum of buckets in the analysis. P-values 
were calculated by univariate ANOVA (a = 0.05) using the normalized values from 
the PCA. As for the qPCR analysis, tests were considered significant when P-values 
lower that 0.05 were obtained. 
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